Hydrodynamic phase-locking of swimming microorganisms 
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Some microorganisms, such as spermatozoa, synchronize their flagella when swimming in close 
proximity. Using a simplified model (two infinite, parallel, two-dimensional waving sheets), we 
show that phase-locking arises from hydrodynamics forces alone, and has its origin in the front- 
back asymmetry of the geometry of their flagellar waveform. The time-evolution of the phase 
difference between co-swimming cells depends only on the nature of this geometrical asymmetry, and 
microorganisms can phase-lock into conformations which minimize or maximize energy dissipation. 



PACS numbers: 47.15.G-, 47.63. Gd, 47.63. mf, 87.17.Jj 

Swimming cells, such as spermatozoa or flagellated 
bacteria, are ubiquitous in nature, yet their dynamics 
in a world without inertia is often counter-intuitive []]]. 
Much insight has been gained in the past about the na- 
ture of the swimming of flagellated microorganisms, from 
the pioneering work of Taylor 0], through to many de- 
tailed reviews on locomotion at the microscale 0, 0, • 
Spermatozoa in particular have received much attention 
in an effort to improve our understanding of the biome- 
chanics of reproductive processes [6[ . 

One particularly puzzling phenomenon observed in 
swimming spermatozoa and other microorganisms is the 
apparent synchronization of the beating of the flagella of 
two or more cells when in close proximity |2|, |2| . This phe- 
nomenon was first modeled by Taylor using infinite two- 
dimensional sinusoidal sheets 0]. Taylor found that the 
most energetically efficient configuration for two swim- 
mers close together was to beat in synchrony. A compu- 
tational model of the same setup showed that at small 
but finite Reynolds numbers, the sheets can achieve sta- 
ble phase-locking at in-phase and opposite-phase config- 
urations 0] , a result which remains valid for finite swim- 
mers with flagella of linearly increasing amplitude Q. 
Computations in two dimensions showed that the flow 
fields of interacting swimmers tend to cluster them to- 
gether into tight synchronized groups 7]. Large arrays 
of cilia (short, closely-packed flagella) also display syn- 
chronization if the internal force-generation mechanism 
generating their beat pattern is load-dependant [lOj, . 

What is still not understood is what constitutes the 
essential physical ingredients to obtain an evolution in 
time to a phase-locked configuration between cells swim- 
ming close to each other. Here we consider a simplified 
model of nearby swimming cells with a prescribed wave- 
form. We show that stable phase-locking can be obtained 
purely passively, due to hydrodynamic interactions. The 
phase-locked state to which the cells evolve is dictated 
solely by the geometry of the flagellar waveforms of the 
cells (specifically, their front-back asymmetry), and not 
by considerations of energy dissipation. 

In the spirit of Taylor Q , we consider a model of co- 



swimming cells consisting of two infinite parallel two- 
dimensional sheets propagating lateral waves of trans- 
verse oscillations with prescribed wavenumber k, fre- 
quency u>, and wave speed c = u)/k] each sheet is thereby 
propelled in the direction opposite to the wave. This ide- 
alized geometrical model, which has been used success- 
fully in the past to study other properties of cell locomo- 
tion [5|, will allow us to clearly elucidate the necessary 
ingredients required for synchronization. We also relate 
it below to experimental observations. 

The shape of the waveform is assumed to be the same 
for both swimmers, and is described by an arbitrary 
function a. The position of the bottom sheet relative 
to an axis about which it is centered vertically, is de- 
noted y± = a(kx — u>t) in its swimming frame. The top 
sheet which is some mean distance h above and paral- 
lel to the bottom sheet moves at a speed Ua relative to 
the bottom sheet. The two sheets have an instantaneous 
phase difference denoted <f> ( see Fig- H); Ua is denned 
to be positive if the upper sheet moves to the right rel- 
ative to the lower sheet; <f> is defined to be positive if 
the upper sheet is left of the lower sheet by <f>. The in- 
stantaneous position of the top sheet is thus given by 
V2 = h + a(kx — u)t — J Q kUA,(t')dt' + <f>), where <j> is the 
phase difference at t = 0. 

The governing hydrodynamics equations for low- 
Reynolds number flow in an incompressible Newtonian 
fluid are the Stokes equations, {Vp = /^V 2 u, V • u = 0}, 
for the velocity field, u = (u,v) and pressure, p. We 
non-dimensionalize using x* — xk, y* — y/h, t* = tu>, 
u* = u/c, U A — Ua/c, v* — v/ec, p* = pe 2 /^,uj, where 
e indicates the ratio of mean separation of swimmers to 
their wavelength, e = hk. Fluid force per unit width 
are nondimensionalized as /* = e///zc and energy dis- 
sipation rate per unit width as E* — e 2 E / fiojch. For 
convenience we introduce the variables h* = y% — y*, 
x* = x* - t*, 4> = 4> — Jo U A (t')dt'. Consequently, 
we have y\ — a*(x*), and y% = 1 + a*(x* + 4>), and 
the phase evolves in time according to <f> = —U A . We 
further assume that the waveform possesses reflectional 
symmetry with respect to the horizontal axis, namely 
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<f> = In-Phase 



(f> = 7T Opposite-Phase 

FIG. 1: (Color online). Our model for flagellar phase-locking: 
Two infinite, parallel, 2D sheets propagate periodic waves at 
a speed c leading to swimming in the opposite direction. The 
top sheet is allowed to move with velocity Ua with respect to 
the bottom sheet and is out of phase by an angle <f>. 



a*(x*+n) = —a*(x* ), in order to focus on cells swimming 
along straight lines (l2|. We now drop the (*) notation 
for convenience, and refer to dimensionless variables. 

As seen experimentally, phase-locking can occur when 
the cells flagella beat close together, therefore an appro- 
priate limit to study is when the mean distance between 
them is much smaller than their wavelength, i.e. e«l. 
In this limit the Stokes equations reduce to the lubrica- 
tion equations, dp/dx — d 2 u/dy 2 and dp/dy — 
We solve these equations in a frame moving with the 
waveform of the bottom sheet [3] • The boundary con- 
ditions are hence u(yi) — —1 and it (2/2) = — 1 + Ua- The 
solution is for u readily obtained as 
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(1) 



Integrating the continuity equation over h yields a rela- 
tion between the gradient of the flow rate between the 
sheets, Q = fjj* udy, and their relative motion as 



9Q _ dy 2 
dx dx 



(2) 



In order to determine the physical conditions for phase- 
locking to occur, we first set Ua — and investigate the 
resultant horizontal hydrodynamic force, f x , acting on 
the upper sheet. In a free-swimming situation, the upper 
sheet would move at a rate such that the viscous drag 
would balance with f x (see below). With Ua = 0, we 
know from Eq. |(5J) that Q is constant and upon integrat- 
ing Eq. flU over h we get Q = — h— (h 3 / 12) (dp/dx). Since 
the system is 2n periodic, we have J Q n (dp / dx)dx = 0, 

which leads to Q = —I2/I3 where Ij = h~^dx. The 
pressure gradient is then obtained to be 



# = 12 

dx 
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(3) 



The force per unit width is determined by integrating 



the stress over the upper sheet, f x 
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FIG. 2: (Color online). Swimmers with reflection symmetries 
by the horizontal and vertical axes cannot phase lock: If a 
relative force exists in [2^,, we obtain the forces in [2b and 
[2J; by reflection symmetries (R planes). Applying kinematic 
reversibility to (KR line) leads to a force in [5J1 which is 
minus the one in [2b, indicating that they both must be zero. 



where n is the unit normal to the sheet into the fluid 
and er = —pi + /x(Vu + Vu T ) is the stress tensor. Using 
integration by parts, the force is given by 
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Using Eqs. HJ and ©, we finally obtain the force 
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<i>) + a(x)} \ dx. (5) 



n ds. 



Physical insight can be gained by inspection of Eq. |5]| . 
When = (in-phase swimming), h is constant, and the 
force is identically zero for all waveforms a(x). Further- 
more, since our waveforms possess reflection symmetries 
about the horizontal axis, the force is also exactly zero 
when <p = 7T (opposite-phase swimming). What is how- 
ever the nature of the hydrodynamic force about the fixed 
points <fi = 0, 7r? 

Using symmetry arguments, we can first demonstrate 
that if the waveforms also possess reflection symmetry 
with respect to the vertical axis (such as a pure sinewave), 
no phase locking is ever possible. This is illustrated in 
Fig. Suppose the force between the swimmers acts 
to stabilize the phase difference (Fig. \2k)- Let us then 
compare the forces on the setup obtained by a reflection 
by the vertical axis (Fig.[2}D), with that obtained first by 
reflection by the horizontal axis and then by kinematic 
reversibility (i.e. change of the direction of the wave) 
(Fig. [2t and d). The force in Fig. [2ji is destabilizing 
while the one in Fig.[5j3 is, for the same setup, stabilizing, 
indicating that both of them must be zero. In particular, 
sine- waves, such as the ones considered in Ref. @, H|, 
cannot phase-lock. Note that this argument holds also 
for finite flagella with a shape invariant upon reflection 
about the vertical axis versus the horizontal axis, such as 
sine- waves with an integer number of wavelengths. 
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If the waveform is not front-back symmetric (i.e. lacks 
reflection symmetry with respect to the vertical axis), 
such as spermatozoa with flagellar waveforms of increas- 
ing amplitude (l5|, the comparison between Fig. and 
Fig. [2Ji cannot be made, and a force can appear. In 
that case, the comparison between Fig. [2^ and Fig. [2t 
shows that this force must be an odd function of the 
phase difference, i.e. f x (—<fi) = —fx(4>)- F° r small vari- 
ations about the fixed points, (p — <f>o + <t>' with 0' <C 1, 
the force given by Eq. ([5|) determines the stability of <j>Q. 
Near the in-phase configuration (</>q = 0) the force is 



fx Q 



-72<p'- 



(6) 



to leading order in <fi' . US denotes the amplitude of the 
waveform, we see that f Xa ~ ±0' 3 £ 4 , with a sign that de- 
pends solely on the wave geometry; the sign + (— ) leads 
to stability (instability) of in-phase swimming. Similarly, 
expanding about opposite-phase swimming (<pQ = 7r), we 
get at leading order 



, 3 [ 2 * (da/dxf 



(l-2a) 4 V(l-2a) J 3 



Ji 



1 dx, (7) 



where J n = ^"(1- 2a) n dx. For small- amplitude waves 
with 5 < 1, Eq. simplifies to 



(8) 



Comparing Eq. © with Eq. ([5]), we see that forces 
near in-phase and opposite-phase configurations have the 
same magnitude, but opposite signs. One of the fixed 
points for phase-locking is therefore stable while the 
other is unstable, in a manner which depends solely on 
the waveform geometry: If the waveform a is such that 
A = f 2 ™ a (da/dx) 3 dx < (> 0) then in-phase swim- 
ming is stable (unstable) while opposite-phase swimming 
is unstable (stable) [lr3 |. 

The rate of energy dissipated in the fluid between the 
sheets per unit width is E = JJ a : Vu dydx, which is 
given over one wavelength by 



E=12 



hh 3 



1 



dx. 



(9) 



For small-amplitude waves near in-phase swimming 
(4>a = 0), we get Eq w I2(j)' 2 L (da/dx) 2 dx, while near 
the opposite-phase configuration (<fi = n), we obtain 
En w 12 J Q 27r (4a 2 — (da/ dx) 2 (j)' 2 ) dx. In-phase swim- 
ming is therefore always the situation where the swim- 
mers have to do the least amount of work, while opposite- 
phase the most work 0, S] 17 1. Comparing this result 
with the forces in Eqs. ((U and ([8|), we see explicitly that 
there is no relationship between viscous dissipation and 
hydrodynamic force, and the swimmers can be forced into 



a stable conformation where the energy dissipation is in 
fact maximum (when A > 0). 

In order to observe the evolution of the phase angle 
towards a phase-locked state, we now allow Ua to be 
nonzero. Evaluating Eq. @ with Eq. |T]) and integrating 
in x we get 



dp 6U A - 12 12U A y 2 + C 
dx 



h 2 



h 3 



(10) 



where C is a constant, found by enforcing that the pres- 
sure is 2tt periodic, C = (6U a -12)(I 2 /I 3 )-12Ua(K/I 3 ), 
where K = J 27r y 2 /h 3 dx. We then obtain Ua by im- 
posing that the swimmers are force-free. Substituting 
Eqs. (JTO]) and ((TJ into Eq. (gj) with f x = 0, we solve for 
<j> = —Ua and get 



dt 



where 



A = 
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1 + a(x + (f>) — a(x) 



dx , 



(11) 

[a(x + 0) + a(x)] 
(12) 



and f x refers to the force in Eq. ((5]). The change in the 
phase angle between the two sheets is therefore propor- 
tional to the force which would be acting between them 
if they were prevented from having any relative motion. 
For small-amplitude waves, we get A -1 « 2ir, so that 
near the fixed points, the phase behaves as 4" ~ ±5 4 <p' 3 
and we have <fi' ~ t -1 / 2 near stable fixed points and 
4>' ~ (t — t)~ x / 2 near unstable points. 

As discussed above, the geometry of the wave is the 
only factor determining the direction in which the rel- 
ative position of the sheets evolve. For illustration, we 
now consider waves in the form of skewed sinusoids. We 
map a sinewave from the intervals [0 : 7r/2] and [7r/2 : n] 
to the intervals [0 : ir/2 + a] and [7r/2 + a : 7r], and define 
a(x + 7r) = —a(x) on the interval [it : 2tt]. Shapes with 
a > (a < 0) have a larger region where the wave am- 
plitude increases (decreases) in the direction of the wave 
propagation (see Fig. [31 upper inset). If a > (a < 0) 
then A < (A > 0) and our analysis predicts that phase- 
locking will occur at the in-phase (opposite-phase) con- 
formation. Experimental observations in Ref. [15j | sug- 
gest the flagellar amplitude of bull spermatozoa is reason- 
ably given by a linearly increasing sine wave; this yields 
front-back asymmetry corresponding to A < 0. We now 
proceed to solve Eq. (fTTj) numerically. 

The dynamics of phase-locking is illustrated in Fig. [3] 
(top) , where we plot the evolution of the phase angle from 
an initial phase <p — tt/2 for three shapes of swimmers 
(upper inset): perfect sinewave (a — 0, black solid line), 
and two skewed sinewaves with a = 37r/10 (blue dotted 
line) and a = — 37r/10 (red dashed line). The perfect 
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FIG. 3: (Color online). Top: Evolution of the phase angle, 
<j>, from an initial value of n/2 for swimmers with three dif- 
ferent waveforms [upper inset]: sinewave (black solid line), 
sinewaves shifted forward and backwards by 37r/10 (blue dot- 
ted and red dashed lines respectively). As predicted by the 
theory, the perfect sinewave shows no phase evolution while 
the skewed sine waves evolve to phase-locked positions (in- 
phase and opposite-phase respectively). The energy dissi- 
pated for the three waveforms is almost identical [lower inset] 
(see text). Bottom: For all initial conditions, the phase an- 
gle converge to a phase-locked state at the only stable fixed 
point, = (same parameters as in top figure, dotted line). 

sinusoidal shape yields no evolution in time, as predicted 
by the analysis. In contrast, the skewed sinewaves evolve 
into the predicted phase-locked positions, in-phase for 
a > and opposite-phase for a < (at the same rate be- 
cause the waveforms are symmetric to each other about 
the vertical axis). The energy dissipated between the 
sheets, Eq. [SJ is displayed in Fig. [3] for the three shapes 
[lower inset]. The dissipation for the two skewed wave- 
forms is identical - as E is invariant under a reversal of 
swimming direction -, is slightly different from the per- 
fect sinusoid (mean square difference < 0.1%), and all 
shapes display the predicted maximum at tt and mini- 
mum at and 2tt. In agreement with our analysis, the 
numerical results demonstrate therefore that swimmers 
with increasing amplitude (a > 0) phase-lock into the 
most energetically favorable conformation, while those 
with decreasing amplitude (a < 0) phase-lock into the 



least energetically favorable conformation. The example 
with a = 37r/10 is further illustrated in Fig. [3] (bottom) 
where we show the evolution of the phase angle from var- 
ious initial phases. In all cases, the configuration evolves 
into a phase-locked state at the only stable fixed point, 
= 0. Further computations (not shown) show that 
increasing the asymmetry of the waveform, or its am- 
plitude, decreases the time scale over which the systems 
evolves into a phase-locked state. 

In summary, in this paper we have used a simplified 
model to show that hydrodynamic forces alone can lead 
to the observed phase-locking between two swimming mi- 
croorganisms if their waveforms are front-back asymmet- 
ric. The nature of the phase-locked state, and the rate 
at which the relative conformation of the two swimmers 
evolve to it, is dictated solely by the geometry of the 
waveforms. In particular, an in-phase conformation may 
be obtained when the swimmers have shapes with in- 
creasing amplitude front to back, as observed for some 
mammalian spermatozoa [To] ]. Other front-back asym- 
metries, such as the presence of a head, would also fur- 
ther contribute to phase-locking Q. 

Discussions with T. Powers and S. Spagnolie and 
funding by the NSF (grants CTS-0624830 and CBET- 
0746285) are gratefully acknowledged. 
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